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Progress by investigators at the Georgia Institute of Technology in the 
development of techniques for passive microwave retrieval of water vapor, 
cloud and precipitation parameters using millimeter- and sub-millimeter 
wavelength channels is reviewed. Channels of particular interest are in the 
tropospheric transmission windows at 90, 166, 220, 340 and 410 GHz and 
centered around the water vapor lines at 183 and 325 GHz. Collectively, 
these channels have potential application in high-resolution mapping (e.g. , 
from geosynchronous orbit) , remote sensing of cloud and precipitation 
parameters, and retrieval of water vapor profiles. 

During the period from January 1, 1993 through June 30, the Millimeter- 
wave Imaging Radiometer (MIR) ocnpleted data flights during a two-month 
long deployment in conjunction with TDGA/CDftKE. Coincident data was 
collected from several other ground-based, airborne, and satellite sensors, 
including the NASA/MSPC AMER, MIT MTS, EKSP SS^/T-2 satellite, collocated 
radiosondes, ground- and aircraft-based radiometers and cloud lidars, 
airborne infrared imagers, solar flux probes and airborne cloud particle 
sampling probes. 

Optimal calibration of the MIR using a Wiener filter to facilitate 
es timatin g the instrument's gain and offset is being implemented for MIR 
data calibration and archiving. Reduction of the clear-air flight data 
shows good general agreement between MIR brightnesses and computed 
brightnesses based on coincident radiosondes, although the discrepancies 
between these brightnesses indicate that the humidity probes on AIR and VIZ 
type radiosondes indicate too moist of a reading in dry regions of the 
atmosphere. 

In order to resolve questions concerning the absolute calibration of 
both the MIR and similar microwave sounders, the investigation of the 
scattering and emission from microwave blackbody calibration loads has 
continued. A new analytical formulation of the coupled wave method for 
lossy period surfaces has been developed, and an extension of the method to 
two-dimensional periodic structures is being developed. A preliminary 
thermal analysis for the wedge-type load is being applied to determine the 
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temperature profiles within such loads. 

Integration of three sulanillimeter-wave channels on the MIR at 325+/- 
1,3, and 8 GHz is proceeding, and should be completed in time to obtain 
data from ER-2 flights planned for September, 1993. 

A study ocmnenoed under this grant is the investigation of the utility 
of the third Stokes parameter for passive remote sensing of both 
ocean wave direction and oriented ice particles in cirrus clouds. Data from 
a polar imetric 92-GHz fixed-beam slant-path radiometer flown on the NASA 
DO-8 during TOGA/OQARE is currently being reduced. Preliminary results from 
constant bank-angle maneuvers show brightness perturbations that are 
correlated with the direction of the aircraft, suggesting that information 
regarding the anisotropy of the ocean surface can be measured using passive 
radiometers. 
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SMART OF ACTIVITIES 


Activities within the period from January 1, 1993 through June 30, 1993 
by Georgia Tech researchers in millimeter and s t itnill 1 m eter wavelength 
tropospheric remote sensing have focused on the optimal calibration and 
archiving of the MIR aircraft data, clear-^air MIR data analysis, operation 
of the MIR aboard the NASA ER-2 during the Tropical Ocean Global 
Atm 3 sphere/Coupled Ocean-Atmosphere Response Experiment (TOGA/OQARE) , and 
preliminary data analysis from TOGA/OQARE. The MIR instrument is a joint 
project between NASA/GSFC and Georgia Tech. In the current configuration, 
the MIR has channels at 89, 150, 183+/-1,3,7 , and 220 GHz. Provisions for 
three addition al channels at 325+/-l,3 and 8 GHz have been made, and a 325- 
GHz receiver is currently being built by the ZAX Millimeter Wave 
Corporation for use in the MIR. Past Georgia Tech contributions to the MIR 
and its related scientific uses have incl u ded basic system design studies, 
performance analyses, circuit and radiometric load design, in-flight 
software, operation during deployment, post-flight data display software, 
and radiative transfer studies of the effects of clouds, water vapor, and 
precipitation on MIR brighli eonos [Gasiewski , 1992] • 

In addition to the MIR-related activities, a 92-GHz polar imetric 
r adiome ter was deployed on the NASA DC-8 during TOGA/CQARE. The purpose of 
this instrument was to investigate the relation between polarized microwave 
brightness tenperatures and both oriented atmospheric ice cloud particles 
and anisotropic ocean surface features. 


l . mir Data Analysis 

All in-flight MIR software for unattended operation aboard the NASA ER- 
2 has been ocnpleted. Post-flight data analysis tools under development 
include a graph ical display program for the PC which incorporates optimal 
calibration filtering. MIR calibration is aoocnplished using hot and cold 
blackbody load views (once per scan) along with appropriately weighted 
temperature values from eight resistive tenperature sensors located on the 
loads. Previous radiometer calibration techniques have used a causal 
infinite iapulse response (IIR) discrete filter to provide estimates of the 
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instr umen t gain and offset values from the noisy single-scan es t i m a t es. 

Recently, software has been written to demonstrate the application of 
an o ptimal non— causal filter (the Wiener filter) to es tim a t e MIR gains and 
offsets [Adelberg et al, 1993]. The filter coefficients are derived frcm 
es timates of the autocorrelation functions far the gain and offsets of the 
instrument. As illustrated using dear-air MIR data (Figure 1) , alcng-track 


brightness variances are reduced by "10-70% using the optimal filter. The 
r eduction s are caused by the removal of calibration noise, bringin g the 
resulting variance doser to the theoretical minimum given by Tgy^/V B?'.' 
Hie improvements obtained using the o pt i mal f ilt e r are most significant in 


the noisier channels (e.g, . 150, 183, 220 GHz); however, even the 90 GHz 
channel, which is relatively stable and quiet, shows inproveraent. 1 The 
software also incorporates ocnpensaticn to correct for aircraft roll. 

MIR data is being analyzed from flights during three ER-2 field 


deployments, sunmarized in Table 1. Applicable data from other instruments, 
including the Massachusetts Institute of Technology's Millimeter-wave 
Ttenperature Saante r (MIS) , the NASA/MSPC Advanced Microwave Precipitation 
Radiometer (AMFR) , the NASA Airborne Ocean Odor Imager (A0d) , the EMSP 
SSH/T-2 radiometric sounder, the NASA/GSPC Raman water vapor lidar, and 
collect*** radiosondes is currently being oonpiled. Several interesting 


r-aco s tudi es have been identified, and are ou t lined as follows: 

(1) Multispectral simulation of the "warm ring” observed at the eyewall 

of tropical cydone Oliver during TOGfl/OCARE. The extent to which this "40- 
km brightness ring is attributable to ocean surface rou^mess or 

low-altitude rain is currently unknown. Resolution of this question is 
e xpected to provide insight into the ability of high-resolution radicmetry 
to differentiate between surface rain and ocea n ic wind, and the ability of 
high-resolution microwave radiometers to monitor the small-scale dynamics 
of tropical storms. Data from the MES and AMER radiometers and the JPL 
AEMAR rain radar on the DO-8 is being sought for use in this study. 

(2) Millimeter-wave radiative transfer simulation of brightnesses 
observed over douds. Investigators at the Kavieng site during TOGA/OQM® 


1 a publication documenting the opti mal filtering technique is being 
p repar ed for submission to the IEEE Transactions on Geoscience and Remote 
Sensing. 
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operated several ground-based instruments for cloud and water vapor 
sensing, i ncludin g ground-based water vapor and cloud radiometers, a cloud 
Ildar, and a cloud oeilcmeter. Data from approximately 3-4 coincident 
overflights of these instruments is being obtained 2 for the purpose of 
conducting ocnpariscns between predicted and MIR-abserved brightnesses over 
clouds. 

(3) dear-air MIR brightness data f rom both airborne and ground-based 
observations is being ocnpared with radiative transfer calculations based 
on nearly coincident radiosondes and Raman lidar data. The results of these 
ocnpariscns (shown in Figure 2 and Table 2) show that the relative humidity 
err o r s present in some types of standard radiosondes produces brightness 
variations (for both satellite- and ground-based observation geometries) 
that are much larger than the attainable accuracy of radio m eters. This 
suggests that relative humidity, when averaged over a altitude range 
ccnparable to the tenperature weighting function width, can ultimately be 
observed more accurately using radiometers than radioso n des. Moreover, the 
ocnpariscns show that the AIR radiosonde provide a consistently positive 
bias in relative humidity, as evidenced by the fact the 183+/-1 GHz channel 
exhibits a large negative discrepancy for the airborne case, but the 150 
GHz channel exhibits a large positive discrepancy for the ground— based 
case. 3 

(4) An interesting case of abrupt upper at m o sph ere drying across a 
convective front has been identified during TOGA/OCARE ER— 2 flight 93-053. 
The f r ra f * ! crossing was also observed by the ER-2 lidar, thus clear-air 
regions can be unambiguously identified, and cloud top levels esti m a t ed. 
The drying ca used an increase in brightness of '*12—14° K in the 183+/— 1 GHz 
channel over a distance of '*'5 Jan? thus this case is expected to be useful 
in studies of the effects of spatial resolution on water vapor retrieval. 


2 Collaborating on the Kavieng overflights is Dr. E. Westwater of the 
NCAA Wave Propagation Laboratory in Boulder, CD. 

3 The results of this investigation being documented for s ub m ission to 
the TFFire Transactions on Geoscience and Remote Sensing. 
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?„ aiifeaafcjgl Tna d Analysis 

A new fo rmul ation of the exiled wave method [Gaylord and Moharam, 
1982] based upon multiport linear network theory has been developed, and 
reflectivity calculations for several wedge abso rb er g eome tries have been 
performed. 4 

In fv pVr to predict the a m o u nt of ther mal emiss ion from a calibration 
target, the temperature distribution over the structure n eeds to be known. 
*ibis is particularly critical near the tips of the absorber, where much of 
the emission occu r s and thermal gradients are large. To this end, a steady 
state thermal analysis of a wedge- type structure has been continued. The 
solution to the heat equation is performed using the ocupled harmonic 
method. Preliminary results suggest that thermal gradients near the 
absorbing (and, hence, emitting) tips of the wedges are significant enough 
to warrant consideration in calibration load analysis. 

3. 325-GHz Receiv er for use cn the MIR 

The implementation of a 325-GHz receiver on the MIR remains a high 
Georgia Tech priority. To achieve an acceptable noise figure (10-12 dB DSB) 
a modif ied design employing a fundamental mixer with DC bias and a 
frequency-tripled 108 -GHz Gunn oscillator is being tested (Figure 3) . 
Fabrication is being carried out by the ZAX Millimeter Wave Corporation of 
San Dimas, GA. A previous design using a double-balanced subharmonic mixer 
was fabricated by ZAX for installation during TOGA/COAKE, but did not 
perform to required specifications . As of this writing, the receiver is 
expected to be available far integration into the MIR near the end of 
August 1993, and operative on the ER-2 for data flights during September- 
October, 1993. 

4. Polarimetric Microwave Radicmetry. 

While the utility of vertically and horizontally polarized brightness 
tenperatures for both surface and atmospheric remote sensing has long been 


4 A paper Axumenting the ocupled wave analysis has been submitted for 
publication in the tfff Transactions on Antennas and Propagation [Gasiewski 
and Jackson, 1993 ; Appendix A] . 
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recognized, the potential use of the third Stokes para me te r 1^=<E V %*> in 
Earth remo te sensing has largely been overlooked. Recent laboratory 
measurements using a fresh-water wave tank [Runkee and Gasiewski, 1993] 
illustrate a significant dependence of % on the propagation direction of 
the water wave at 92 GHz. Moreover, the angular dependence complements that 
exhibited by Ty and Tfc, suggesting that passive remote seising of water 
surface wave direction can be facilitated by polarimetric microwave 
radicmetry. As suggested by the work of Wentz [1992], polar imetric 
radicmetry could have potential application in measuring ocean surface 
winds. Other investigations by Dzura et al [1992] suggest that polar imetric 
radicmetry will also be useful far observing nonlinear ocean wave 
interactions. In addition to striated water surfaces, it is hypothesized 
that polar imetric microwave signatures in T v , Tfo, and T\j will be produced 
by oriented ice particles, for example, in electrified cirrus anvils. 
Indeed, microwave depolarization signatures in space-to-ground 
oonnunications links have been associ with lightning disch arges [Oooc 
and Arnold, 1979] . 

To develop a theoretical model for the polar imetric surface emission 
process, we have implemented a geometrical optics model for the upwelling 
Stokes' vector. The model incorporates double scattering for steep 
observation angles. Using this model we have corroborated the 
experimental data for the case of a 65° observation angle (Figures 4a and 
b) . The encouraging resul ts of the laboratory wave tank experi m e n t prompted 
an investigation into the practical application of polarimetric radicmetry. 
Necessary electronic and hardware modifications to operate the NASA/GSPC 
9 2 -GHz polarimetric radiometer on the NASA DO-8 aircraft during TOGA/OQARE 
were performed. The radiometer, developed under a previous NASA grant (NAG 
8-829), permits precise measurement of the first three Stokes' pa r a me te r s. 
Accurate pre- and post flight calibration data was collected using the 
pr>i ar-i ypd calibration load described by Gasiewski and Runkee [1993]. A 
total of twenty DO-8 flights were made during TOGA/CQARE, including seven 
lew— altitude (~2-4 km) constant bank-angle maneuvers designed to provide 
constant observation angle views of the ocean surface over a range of 
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azimuthal angles. 5 

pest-mission calibration of some of the 92-Qiz polar imetric data is 
ccnplete. A plot of the constant bank-angle data far flight 93-01-13 shews 
residual peak-to-peak brightness variations with azimuthal angle for Ty and 
% of amplitude 2.5-6 K (Figures 5a and b). The beam observation angle 
relative to nadir is 65°. The solid and dashed curves indicate oos(2J<) and 
sin (2^) variations, (respectively) , where (p is the az imut hal angle of 
observation with respect to the apparent wave direction. This data, along 
with the laboratory wave tank measur ements corroborate the hypothesis that 
anisotropic ocean surface features can be observed using microwave 
radiometers. 


5 This is believed to be the first fully polarimetric radio m e t e r to be 
operated in Earth remote sensing experiments on a US research aircraft. A 
predecessor 37-GHz instrument was flown during experiments conducted by the 
Soviets during the late 1980's [Dzura etal, 1992]. 
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OCNCUSIGNS AND FUNS PCR FUIUFE WCRK 


To date, Georgia Tech has provided assista nc e to NASA/GSPC in the 
mechanical and electrical design of the MIR, and has been primarily 
responsible far MIR flight and analysis software. In addition, the 
investigations conducted ever the past six months (1) confirm the utility 
of Wiener filtering in radiometric calibration, (2) illustrate significant 
humidity biases associated with certain types of radi o so n des, particularly 
the AIR sondes, (3) have provided several case studies using MIR during 
TOGA/OQARE, and (4) suggest the utility of fully polarimetric microwave 
radiemetry far identifying ocean surface anisotropies, with the potential 
application to inproving measurements of ocean surface wind direction. 

l. MIR Observations and Data Analysis 

At this time, our primary interests are the integration of the 325-GHz 
receiver into the MIR, EHF and SMW radiative transfer stu d ies using MIR 
data, and the development of EHF and SMW meteorological retrieval 
algorithms. The EHF case studies will rely mostly on MIR data observed 
during TOGA/CQARE and previous MIR flights, although new data, particularly 
at 325 GHz, is critical far both EHF and SMW analyses. 6 At least one 
Georgia Tech investigator will be present to assist in the MTR deployment 
at Wallops Island during the September 1993 flights. 

The highest Georgia Tech priority is integrating the 325-GHz receiver 
into the MIR for data flights during this deployment. Collaboration with 
ZAX will be made as necessary to help insure that this receiver is finished 
and available for the September flights. Alternate designs far the 325-GHz 
receiver are currently being studied in the event that the ZAX design 
cannot provide acceptable noise figures for airborne radiometric studies. 

All MIR data is currently stored in raw format on high-density 8 -ram 
tapes . To facilitate meteorological data analysis and provide a practical 
means of disseminating MIR data to collaborating investigators, the display 


6 Due to delays in fabricating the 325-GHz receiver, no SMW data is 
yet available. A request to extend this grant through 1994 has been 
submitted so that such data can be obtained and analyzed. 
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so ftware is being modified far disk-based storage. Ibis software, which 
will be available for use cn PC's, will also incorporate a variety of 
siiple user— interactive features for KIk data analysis. data is 

being archived cn 35-nm slide s far gr^hical storage and 
dissemination. 

IWo of the MIR case s tudi es identified from TOGA/OQAKE are particularly 
interesting for radiative transfer experiments, namely, the overflights of 
cyclone Oliver and the Kavieng ground observation site. To this end, it is 
planned to simulate ocean surface and atmospheric conditions within the eye 
of cyclone Oliver to determine the consistency of wideband radiative 
transfer and the observed brightness data. The radiative transfer 

Tnrrfoi will consider the effects of increasing humidity, precipitation, and 
ocean roughness near the eyewall to determine the relative contributions of 
these components to the warm ring. This expe rim ent will require sane data 
from the NASA/MSPC AMER and JPL AFMAR, both of which are available and are 
being requested. 

The Kavieng overflights will be used to provide data for sinulating the 
effects of clouds on upwelling EHF brightness temperatures. Ground-based 
radiometer data will be used to determine total water vapor and cloud 
content, while ground-based lidar will be used to determine cloud bottom 
altitudes. 7 Radiosondes will be further used to constrain the vertical 
distribution of water vapor. Of interest is a comparison of computed and 
observed upwelling brightness temperature, and in particular, the impact 
that clouds have on these temperatures. 

In addition to validating radiative transfer models for TOGA/ttlARE case 
studies, high-resolution retrievals of water vapor and cloud water content 
from EHF and SMW observations will be investigated, ttie premising method 
is the non-linear statistical iterative method used by Kuo [1988]. An 
extension of this technique to incl u d e cloud water content, base altitude, 
r i n ii thickness and surface emissivity as variable parameters is currently 
being studied. Another promising method is based on a multidimensional 
no nlinear map implemented using an artificial neural net. Here, the use of 


7 Airborne l idar data from the ER-2 CIS is not available for the 
Kavieng overflights. 
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the backprcpagaticn algorithm to train a net ocnsisting of two-layers (of 
about 10 neurons each) on simulated radiometric data is currently being 
studied, and will be continued. 

3. Radi meter Calibration 

In order to better characterize the RF response and absolute accuracy 
of the MIR, we plan to conduct several tests on the MIR, i nclu d in g: (1) RF 
passband response measurements using a pi*gma discharge noise source and IF 
spectrum analyzer, ( 2 ) local pprH 1 1 wt-m* interference and reflection 
measurements using a stepped reflecting plate, and (3) calibration-load 
foam reflection and transmission measurements. These simple measur ements 
will provide answers to questions concerning the calibration of the MIR and 
the use of the 183 and 325-Giz data in radiative transfer interccnparisons. 

Accurate absolute calibration of the MIR requires that the total 
reflectivity of the hot and cold loads be less than 1%, and known to better 
than 0.1%. Manufacturer's specifications typically provide only the 
ppyaiiar oonponent of the reflectivity, which is thought to be 
substantially less than the total reflectivity. In order to refine the MIR 
calibration, we plan to extend the study of the electromagnetic 
characteristics of wedge-type blackbody loads to the more desirable 
pyramidal leads. 

Our approach is to develop numerical models for one- and two- 
dimensicnally periodic lossy gratings using the coupled wave method. We 
currently have software based on the coupled wave method to predict the 
reflectivity of one— dimensionally periodic loads of arbitrary d i el e ctric 
profile. Extension of the coupled wave model to tvo-dimensionally periodic 
surfaces appears feasible, and will be attempted. In addi tion to the 
electromagnetic analysis, work on the steady state thermal analysis for 
one— and two-dimensionally periodic calibration loads will be continued. 
Although, the radiometric calibration load reflectivity analysis is of 
in understanding precision radiometer calibration, it is of 
tertiary relative to the EHF and 9MW data collection and 

analysis. 
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3. PolariingfaHc RadirmPfhv 

Further investigation of polarimetric radicmetry under this gran t will 
consist of analysis of 92-01 z polarime tr ic data observed during both 
TOGVOCARE and the laboratory wave-tank and striated soil and sand 
experiments. The limitations of the geometrical optics model in explaining 
the wave tank and ocean surface data will be of greatest interest; the 
validity of this m od e l at 92-GHz will provide insight into the 
applicability of geometrical optics-based models for higher microwave 
frequencies. Information on the polarizing properties of surfaces will be 
useful in remote sensing of ocean characteristics and vegetation as well as 
underst an di n g the effects of surface emission on passive atmospheric 
sounding. 

To r educe the oonplexity of calibrating a polarimetric radiometer, the 
design of a digital cross-correlator will be continued. The cross 
correlator will be a precursor to the one p roposed to be used on the 
NASA/M5FC Advanced Microwave Precipitation Radiometer (AMER) [Gasiewski and 
KUnkae, 1992b] . The bandwidth of the cross-correlator will be approximately 
500 MHz, making the device useful for wideband radiometric channel studies. 
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Table 1 


HER integration and data flights: 1992-93. 

Sortie # Date Time Instruments 

(UTC) & Remarks 


92-087 

5/11/92 

1900-2200 

92-089 

5/14/92 

2315-0515 

92-090 

5/15/92 

2320-0500 

92-130 

7/23/92 

2100-2315 

92-131 

7/29/92 

0700-1330 

92-132 

7/30/92 

0700-1330 

92-134 

8/2/92 

0700-1030 

92-135 

8/3/92 

0700-0615 

92-140 

8/6/92 

0700-1330 

93-053 

1/12/93 

2130-0430 

93-054 

1/17/93 

2300-0710 

93-055 

1/19/93 

0130-0922 

93-056 

1/25/93 

2315-0700 

93-057 

1/30/93 


93-058 

1/31/93 

2200-0600 

93-060 

2/5/93 

1430-2050 

93-061 

2/7/93 

1555-2115 

93-062 

2/9/93 

1815-0025 

93-063 

2/10/93 

1430-2225 

93-064 

2/19/93 


93-065 

2/20/93 

1900-0335 

93-066 

2/22/93 

1900-0205 

93-067 

2/24/93 

2000-0315 


MIR, MIS (u) ,AOCI 
MIR, MIS (d) ,AOCI * 

MIR, MIS (u) ,ADCI * 

MIS(d) 

MES(d) +* 

MTS(d) +* 

MTS(d) + 

MIS(u) + 

MIS(u) +* 

Radiation (93-01-06) 
Convection (93-01-07) 
Convection (93-01-08) 
Radiatiorv/Kavieng (93-01-09) 
Pilot proficiency 
Radiation (93-01-10) 
Convection, Oliver overflight 
(93-01-11) * 
cyclone Oliver, MES(d) 

Oliver overflight (93-01-13) 
Kavieng, HES(d) (93-01-14) * 
Pilot proficiency, MIS(u) 
Convection (93-01-16) 
Convection (93-01-17) 
Radiation (93-01-08) 


* SSJVT-2 satellite underpass. 

+ Ground-based Raman HoOy lidar overflight, 
(xx-xx-xx) indicate coincident DC-8 flights. 

Note: "u" or "d" indicate vqp-looking or dcwn- 
looking, respectively. 
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Radiosonde 



Frequency (GHz) 



Instrument Type 

89 

150 

183±1 

183±3 

183±7 

220 

AIR 

Mean (K) 

3.73 

2.78 

-11.49 

-6.92 

-3.27 

0.00 

Std. Dev. (K) 

2.82 

3.45 

3.61 

3.49 

2.27 

1.36 

# samples 

18 

18 

18 

18 

18 

18 

VIZ 

Mean (K) 

3.93 

1.85 

-6.93 

-2.63 

-0.67 

-0.03 

Std. Dev. (K) 

2.84 

1.02 

1.59 

1.57 

0.94 

0.61 

# samples 

9 

9 

9 

9 

9 

9 

Viasala 

Mean (K) 

-1.98 

-1.04 

-2.70 

-0.20 

0.79 

0.53 

Std. Dev. (K) 

2.07 

1.80 

3.06 

1.88 

1.34 

0.76 

# samples 

12 

12 

12 

12 

12 

12 

Lidar (w/AIR) 

Mean (K) 

0.22 

0.94 

-11.65 

-5.05 

-1.83 

0.15 

Std. Dev. (K) 

2.30 

3.16 

3.20 

2.47 

1.89 

1.94 

# samples 

125 

125 

125 

125 

125 

125 

Lidar (w/VaiBala) 

Mean (K) 

0.15 

0.52 

-6.48 

-2.94 

-0.63 

0.06 

Std. Dev. (K) 

1.80 

1.83 

2.16 

1.43 

1.08 

0.86 

# samples 

69 

69 

69 

69 

69 

69 


Table 2. Comparison of nadiral airborne ocnputed and observed 
brightness for KIR channels rising AIR, VIZ and Vaisala-type 
radiosondes. The data are from the Wallops Islands aircraft 
flights (see Table 1) . 
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Variance Variance . Variance 


89 GHz 150 GHz 



Pixel Pixel 


Figure 1: Briffitness variance (in K 2 ) vs. scan angle (i.e., pixel 
number) for MIR data calibrated using an IIR filter (solid 
curves) and an optimal Wiener filter. Variance reductions of ‘10- 
70% are obtained vising the optimal filter. 







Record Number 


Figure 2 . Ocnpariscn of ground-based oaqputed and observed 
brightness for MIR channels using AIR and Vaisala-type 
radiosondes. The data is from Wallops Island on 7/30/92. 
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Reduced Height 
RF Waveguide 


LO Backshort 



Whisker 

Contacted 

iode 



Filter to Block RF 


Quartz Substrate in channel 


DC-Bias Input 


H — IF OUTPUT 


Low Pass Filter 
to Block LO and Pass IF 


XSNNSs LO Waveguide 


Figure 3. 325-Ghz wHver configuration, using waveguide Schottky 
ri-i nfte nixing with a fundamental local oscillator. 
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Brightness Temperoture (K) Brightness Temperature (K) 





- <Th> (K) Tv - <Tv> (K) 


TOGA/COARE 93-01-13 23:58:40 - 00:02:23 Z 



Beam Azimuth Angle (deg, relative to obs. wave direction) 


(a) 



(b) 


Figure 5. Vertical (a) and horizontal (b) 92-GHz brightness 
perturbations observed during constant bank-angle turns over 
ocean during 70GA/OCAKE. Hie beam angle is “65° free nadir. 
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Abstract 


The scattered, transmitted and internal fields of a lossy periodic absorber il- 
luminated by a plane wave are analyzed using coupled wave theory. The solution is 
facilitated by a cascaded linear-multiport formulation in which each layer of the peri- 
odic structure is considered to be a linear multiport network. A method for determin- 
ing the required number of harmonics for acceptable convergence of the coupled-wave 
field expansions is developed. 

The coupled wave reflectivity calculations are corroborated by measured data 
obtained from millimeter-wave bistatic scattering measurements of wedge-corrugated 
radiometer calibration targets. The angular response of the scattering range antennas 
is remove by deconvolution using the CLEAN algorithm. The coupled wave reflectiv- 
ities are also shown to be consistent with calculations based on an integral equation 
method. 

Using the coupled wave method, design data for wedge-corrugated iron-epoxy 
calibration targets is presented. The calculations show that total reflectivities of less 
than -30 dB over a decade in bandwidth are practical. 
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1 Introduction 


In the design of blackbody targets for microwave radiometer calibration, the scat- 
tering and emission characteristics of periodic lossy dielectric structures are of great 
interest. For precise absolute calibration, it is required that the radiated brightness 
temperature of the surface be known to within ±0.1 K or better. This, in turn, re- 
quires that the emitted radiation intensity be known to within ±0.03 % or better. 
Given the uncertainty in the brightness of typical background radiation fields, this 
also requires that the total fraction of radiation reflected from the target, or total 
reflectivity, be extremely small (often less than -30 dB). Since the total reflectivity 
accounts for scattered radiation that can be incident from any angle, this quantity is 
larger than the commonly-measured retroreflectivity, or backseat tering cross-section. 

In determining both the emitted radiation intensity and the total reflectivity, 
it is essential to know the full solution for the electromagnetic fields scattered from, 
transmitted through, and internal to a periodic structure illuminated by an incident 
plane wave. From this solution, the emitted intensity can be obtained using detailed 
balancing, along with knowledge of the structure’s temperature distribution. Like- 
wise, the total reflectivity can be obtained from the full solution by computing the 
bistatic scattering coefficient of the structure. The full solution is also useful in the 
design of anechoic material for RF absorption and frequency selective surfaces used 
in antenna components, microwave filters, optical dielectric gratings and acoustics 
filters. 


We restrict ourselves in this paper to rigorous numerical solutions to the scat- 
tering problem, that is, those involving no analytical approximations to Maxwell’s 
equations and the electromagnetic boundary conditions. These solutions are exact 
in the limit of infinite computational precision, memory and time. In practice, error 
in a rigorous numerical solution occurs only as a result of finite arithmetic preci- 
sion and limited array sizes. Thus the solution is applicable to surfaces of arbitrary 
electrical scale. Indeed, electrical sizes in between the regions of applicability of geo- 
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metrical scattering (at short wavelengths, [1]) and the extended boundary condition 
method (at long wavelengths, [2]) are of interest in calibration target design. For 
such intermediate-sized surfaces the Rayleigh criterion for smoothness [3, 4, 5] is not 
generally satisfied. As a consequence, the surface exhibits a complex spectrum of 
both upward and downward travelling waves within the troughs. 

Several rigorous numerical methods for periodic structure analysis have been 
developed. These fall under the general headings of integral equation methods using 
periodic Green’s functions [6, 7, 8], finite difference methods [9], finite element meth- 
ods [10], coupled wave methods [11] and variations or combinations of these [12]. In 
previous studies of microwave absorbers [1, 13, 8], only the retroreflectivity of the 
periodic surface had been considered. In this study, we consider both the total reflec- 
tivity (as in [6]) and the internal and transmitted fields. Due to its applicability to 
to general periodic structures which might contain voids, dielectric layers, and other 
inhomogeneities, the coupled wave method is preferred. 

Rigorous coupled wave solutions to the TE case have been formulated for 
corrugated periodic surfaces [14] and to the TM case for non -corrugated dielectric 
gratings [15]. However, both of these formulations have relied upon solutions to inho- 
mogeneous second-order vector wave equations. A more direct solution, particularly 
for the general incidence case, has been based upon the first-order curl equations, (c.f., 
Eqs. 14 and 15), but formulated only for non-corrugated dielectric gratings [16, 17]. 
This paper considers both TE and TM formulations for corrugated periodic surfaces 
based on the curl equations [18]. The formulation is symmetric in ~E and 77 , and 
its solution is based on a simple analogy with a cascade of linear multiport networks 
(e.g., [19, 20]). Using the cascaded network formulation and employing generalized 
scattering matrix theory the full field solution can be efficiently computed. 

The coupled wave theory is applied to the analysis and design of wedge- 
corrugated microwave absorbers. This one-dimensionally periodic structure is com- 
monly used in microwave radiometry as a calibration target. Both the internal fields 
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and total reflectivities are shown. The calculations are verified against millimeter- 
wave laboratory bistatic scattering measurements. Due to the relatively wide angular 
response of the scattering range antennas, deconvolution of the measured data us- 
ing the CLEAN algorithm [21, 22] is performed. The coupled- wave reflectivities also 
compare favorably with calculations performed using an integral equation method 
based on a periodic Green’s function. Since it is often desireable to use only a single 
wideband target to calibrate a set of channels spanning over a decade (or more) in 
bandwidth, computed data for the design of wideband corrugated absorber surfaces 
based on the coupled wave calculations is presented. 


2 Periodic Surface Scattering 


The periodic surface scattering geometry for an obliquely incident, arbitrarily polar- 
ized uniform plane wave is illustrated in Fig. 1. The incident wave is 

% = ’E 0 e~ il '- 7 (1) 

where the wave vector k t = k x x + k y y + k t z. The surface grating vector is 7T = Kx, 
where K = 2tt/A and A is the period of the surface. The two homogeneous half- 
spaces (regions 1 and 3) are characterized by complex relative material parameters 
and e 3 ,/t 3 , where e = d - je" and ft = /i' - jft". The periodic inhomogeneous 
region exists from z = 0 down to z = —d. 

The time harmonic Maxwell’s equations for source-free media are: 

Vx5 = — ju>n 77 V x 77 = jut ~E 

V-iE = 0 V • fiE = 0 (2) 

where both e(z,z) and ft(z,z) are periodic functions of x. Since the exciting field 
is also periodic in x, we can use the Bloch-Floquet theorem [23] to express the field 
solutions in the form: 

F(x + = e~ ,ksA F(x,y,z) (3) 


4 



where ~F = ~E or 77. A complete solution for the fields in regions 1 and 3 can be 
written as a series of Floquet harmonics: 

E,(x,y, 2 ) = E i+ £ ft.e-**'-* 4 "** 1 -' «) 

n=— oo 

E>(*,V,2) = £ (5) 

n=— oo 

where 7£„, 7 n are unknown amplitudes. Due to phase matching across region 2, 
identical space harmonics of x and y are used in both regions 1 and 3. The Floquet 
condition (Eq. 3) requires that 

k xn — k x + n/f. (6) 

Since the field solutions must satisfy Maxwell’s equations, we also have 

= -Jk] - Wn + «=1.3 ( ? ) 

where fcj = (w/cj^el/ir and * 3 = The associated magnetic fields can be 

obtained from Eqs. 4 and 5 using Ampere’s law. 

Each Floquet harmonic represents a plane wave mode, the summation of which 
form a complete expansion for the scattered and transmitted fields. Only a finite 
number of these are propagating uniform plane waves; a countably infinite number of 
them are nonuniform evanescent plane waves. For the general case in which the plane 
of incidence does not contain the grating vector (i.e., fa ^ 0, *) the wave vectors for 
the scattered propagating plane waves lie along a cone. In region 1, the polar and 
azimuthal angles for the propagating scattered wave are: 

0, = tan-Mfe + k 7 y )/k lzn ) (8) 

fa = tan _1 [A: v /A: rn ] (9) 

If either region 1 or 3 is lossless, the associated value of k qtn is either positive real 
(k qxn = k' qzn ) for propagating waves or negative imaginary (k qzn = -jk' qin ) for evanes- 
cant waves. If either region 1 or 3 is lossy, the values of k q and k qtn are complex, but 
k xn and k v remain purely real. 
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For the special case when 4> t = 0 or v, the implications of the Floquet condition 
and phase matching are illustrated in Fig. 2. Propagating harmonics in regions 1 and 3 
have associated wave vectors whose tips lie on two circles with radii equal to either 
or k$ (respectively). For the evanescent harmonics, the imaginary component of each 
k 9xn is given by one of two hyperbolic curves, defined by — k^ zn = k*, q = 1,3. 


2.1 Coupled Wave Analysis 


In the coupled wave analysis for periodic surfaces [16, 17], region 2 is discretized 
into L inhomogeneous layers of equal thickness Az (Fig. 3). Within the layer, the 
complex relative permittivity and permeability can be represented approximately by 
known Fourier series: 


e(x, z/ — Az < z < zi) 

= €/(*) = 

E a n/e jKnr 

(10) 



oo 


/j(x, z\ — Az < z < z/) 

= n,{x) = 

E h ^ Knx 

(11) 


n=— oo 


where z\ is the top of the layer, l = 1,2, ...L. The electric and magnetic fields 
within each layer can similarly be represented by expansions in space-harmonic fields: 

^2 i(x,y,z) = E {E xln (z)x + E vln (z)y + E zln {z)z) e-* k '"*+ k ^ (12) 

n =“00 

H 2l {x,y,z) = E {l/ ll nWi + ^(:)y + ^»Wi}e- i( ‘‘" ,+ *’ ,) (13) 

n=— oo 

where the E x t n , E v t n ,E z i n and H x i n ,H v i n , H x i n are the unknown amplitudes of the 
space harmonic. 

Substitution of Eqs. 10-13 into Faraday’s and Ampere’s laws yields a set of 
differential equations in which all space harmonics of the electric field are coupled 
with all those of the magnetic field by virtue of the periodicity of c and \i\ 

V x ~Eu = -jufi 0 m(x)H 2 i (14) 

Vxl 2 i = jut 0 ti(x)E 2 i (15) 
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The above equations can be harmonically balanced [24] to yield an infinite set of first- 
order coupled differential equations for the space harmonic amplitudes (the coupled 
wave equations): 



\E xl ] 


rz-i 

d 

dz 

Hyi 

H xl 

= A 

Eyl 

H xl 


Hy, J 


. Hyl . 


(16) 


where 

(17) 

and ~E x i,~H z i are determined directly from ~E x i,~E v ; and H x i,TIyi, respectively. The 
submatrices and Ai\ can be computed in terms of known parameters such as 
u},K,^i,a nh and b nl . 


0 

Au 

A 21 

0 


For the general incidence-angle case (<f> ^ 0, ir), all of vector components of 
the the space harmonic amplitudes are coupled through the matrix A. However, for 
the special case of <f> = 0,7r, the space harmonics produced by TE-polarized (E 0 = 
E 0 y ) incident waves are decoupled from those produced by TM-polarized (H 0 = 
// 0 y)incident waves. In this case, a complete solution can be found by consideration 
of the TE and TM polarizations separately. 

Consider a TM-polarized incident field with 7f t - = H 0 ye~^ k * x+k,z \ for which 
only H y i,E x i and E z i will be nonzero. Upon harmonic balancing Eqs. 14 and 15 
become: 

d 00 

~ -J k *nE z l n (z) - jT) 0 k 0 (18) 

m=-oo 

-j-Hyln(z) = ~j— J2 E *hn( 2 ) a m-n,t (19) 

dz Vo m= -oc 

E tin (z) = -T } 0 f; tj- Hy lm (z)a m - ntl (20) 

m=-oo 

for n = — oo ... oo, where the a n / are the Fourier coefficients for l/c/(x), k 0 = 

and T} 0 = Since the E z i are linearly related to H y i , only two state variable 

functions ( E x and H v ) are required to completely specify the fields, and the coupled 
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wave equations reduce to the following first-order constant-coefficient linear system: 

( 21 ) 


d 

[Erfl 

[ 0 

Wl 

*-* 

*0 

dz 

Tfy, 

An 

nrj 


T xl 

Uyl 


where 


{^u} nm = jVoko O n -m,l ~ ^n-n,l| 


{^ 21 } = -j—a n -m,i 

y. Jnm rj o 


( 22 ) 

(23) 


In practice, the above system is truncated at the N^ 1 harmonic. Thus, the sub- 
matrices A 12 and A 2 \ are each of size (2 N + 1) x (2N -f 1). 


A solution to Eq. 21 can be obtained by eigenanalysis, in which the space- 
harmonic amplitudes are expressed as a linear combination of eigenfunctions: 

Uyl 


= W Z(z) V 


(24) 


where 


f(z) = 


r -*i« 


0 


0 


0 


0 


(25) 


0 • • • 0 e Aa t aw + 1 )* 

In the above, the A’s are the eigenvalues of A for the fa layer, the columns of W 
are the associated eigenvectors, and V is a vector of 2(2 N -f 1) coefficients, to be 
determined by boundary conditions on the tangential ~E and U fields at the top and 
bottom of each layer. 


2.2 Linear Multiport Solution 


In the transmission matrix solution, each layer of the periodic structure is treated 
as an isolated linear multiport network (Fig. 4), with the layer’s transmission matrix 
determined using coupled wave analysis. The transmission matrix for the entire 
periodic surface is subsequently obtained by multiplying, in sequence, the individual 
transmission matrices to form a multilayer transmission matrix. 
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Consider the incoming and outgoing plane wave modes for a single grating 
layer where both evanescent and propagating modes are included. For the TM case, 


the electric fields are 

represented as: 


I^(x,z) 

_ ^ £ + ^ rn £ j e ~i( k * 

n=-oo l ko ko J 

(26) 

N 
- 1 

N 

= £ '" XT ’ _ ‘ST*/ 

(27) 

l£(x,z) 

= £ E tn j - ^z} e- j( **" x+W) 

(28) 


n=-oo l*® k ° > 




(29) 


where the superscripts refer to incoming (+) or outgoing ( ) waves, and 
k zn = yjkl — kl n . The associated magnetic fields are: 


Ttf(x,z) = £ — ^ y (30 ) 

n=-oo 

THx,z) = £ (31) 

n=-oo Vo 

= £ Sij )e- j ^” x+k ‘ nZ) (32) 

n=— oo Vo 

= t — ( 33 ) 

n=-oo Vo 


The transmission matrix T t for the 1^ layer is defined as: 



Tp 

£ 3 " ' 


' T„ 

r 13 ‘ 


£3" 

l^r J 

= Tj 

[#] 


.5*31 

T33 _ 


J 


where , ~E^, and are vectors containing the plane wave amplitude coef- 
ficients. The four sub-matrices of T t are generalized transmission (matrix) coeffi- 
cients [19]. They play the same role as the familiar scalar transmission coefficients in 
cascaded microwave circuit analysis. 

By harmonic balancing, the tangential electric and magnetic field boundary 
conditions for the top (z = Az) and bottom (z = 0) of the grating can be expressed 
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in matrix form as: 


~%TMo (E± + ~Ei) 


<-*? + 17) 

ZrMoCEt + 


= i, 0 W?(Az)C 


= V.WV 


(35) 


(36) 


on-K) 

where ~Ztm 0 is the TM wave impedance matrix for free space, with diagonal elements 
k 0zn rj 0 ! k 0 and zeros elsewhere. The generalized transmission matrix coefficients are 
subsequently determined by solving Eqs. 35 and 36 under the conditions that = 0 
(for 7 U and jT 31 ), and, alternately, EJ = 0 (for T 13 and T^). These steps lead to: 


= x 


Z~tmA p ee{&z)Z T Mo ~ - Phe(Az)Z T m 0 + Thh(Az)^ 


T 13 = r 


7* = r 


Tz 1 = ^ ^Zj Mo {Tee{Az)Z T m 0 - 7eh(Az)) + T H e{Az)Z T Mo - 7«tf(Az)] 
'%tmo(Pee{&z)Ztmo + ~7eh{Az) ) ~ T*He{& z )~%TMo — 7 * W //( Az)j 

(Az)] 

(37) 


Z T m 0 (Pee(^ z )^TMo + 7 E h{Az)) + 7 H e(Az)'Z TMo + ~?hh{ 


where the P {j are the four sub-matrices of the propagator matrix T(A z) (e.g., see [25]): 


P(Az) = 


Pee(Az) 

T E h{Az) 

_ Phe(Az) 

Phh{Az) 


= We{Az)W 


-1 


(38) 


To find the overall transmission matrix for a multilayer structure, the trans- 
mission matrices for the individual layers are cascaded [19, 20]: 

f = T 1 T 2 T 3 ...Tl-i¥l (39) 

Note that 7 is computed under the assumption of i? 0 as the normalizing characteristic 
impedance. To determine the transmission matrix for arbitrary bounding media, we 
must re-normalize to the impedances of the upper and lower half spaces. This is 
accomplished by first computing the impedance matrix ’Z of the structure, defined 
implicitly by: 


Z TMo /v 0 (Et + ) 

_ 7 

(-Ei + ^1 )ho 

ZtMo/Vo (£3 + ) . 

— Lt 

. (Et - "EDh. . 


( 40 ) 
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In terms of 7 : 

I = 


I ZtMo 

LL_ 1 


1 

Zj 

'13 

Mo J 


Z 33 



^33 + 7i3 

7*+7,il f" 7^33 7*13 

— Tn 

7 

7 J [ 7 

-7 


-l 


(41) 


From the impedance matrix, the scattering matrix 3? of the multilayer structure, 
defined by: _ _ 

#1 


pri 



S 13 ] 

. e ; . 


^31 

^33 . 


(42) 

can be computed from Eq. 40 for arbitrary imbedding impedances rji and 77 3: 

3 = (5 + 7)-* (I - 7) (43) 

where: 


Z = 


—Z TM1 Zn 

~^TMi Zl3 (*?l/*?3) 


ZjM3 Z33 


(44) 


and ~ZjMq is the TM wave impedance matrix with diagonal elements k qzn T) q /k q . The 
scattering matrix consists of four sub-matrices 3>„. These are the generalized scat- 
tering (matrix) coefficients. 


With respect to Eqs. 4 and 5, the amplitudes R n and 7 n for an incident TM 
plane wave are: 


7L 2 

TLo 

7* 2 

_ 7J 

0 

0 

y\H 0 

0 

0 

• 


. 

7- 2 


0 

T.i 


0 

7_o 


0 

7i 


0 

7 2 


0 

• J 


m 9 m 


(45) 
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To explicity compute the fields within an arbitrary layer, the coefficients for 
that layer can be found by stepping downward from the top boundary. The recursion 
is: 


= e (+ \ ( Az /+ i ) W , + \ W i V t 


(46) 


where for the top layer: 


v, = 


TZtmi (ffi + ^7) 1 

(-K + XT) . 


(47) 


and ~E is computed via Eq. 42. Equivalently, stepping can be performed from the 
bottom layer upward. To minimize the accumulated effects of roundoff error, it is 
best to step from both boundaries (z = 0 and z = — d) inward to z = —d/2. 


The above linear multiport formulation has advantages over previous coupled 
wave formulations [14, 15, 16, 17] in both computational complexity and physical 
interpretation. By cascading the transmission matrices for L grating layers the num- 
ber of calculations in the above follows ~ LN 3 . Formulations for which all bound- 
ary conditions between layers are prescribed simultaneously result in complexities of 
~ (LN) 3 . An exception is [17], in which the complexity was reduced by exploiting the 
sparseness of the boundary condition matrix. However, these formulations obscure 
the physical nature of the fields. In the linear multiport formulation the electromag- 
netic fields above and below each layer are simply interpreted as superpositions of 
reflected and transmitted Floquet harmonics. 


Although we have only discussed the TM solution, the TE solution is readily 
found using duality [24], where ~Etm ~Hte> ~Ete, and £«-►/*• Note that 

in the case of a homogeneous layer, the submatrices in A (Eq. 21) are diagonal, thus 
the eigenvectors matrix W is also diagonal, and the eigenvalues are the k gxn . Both 
the propagator and transmission matrices are diagonal, indicating coupling of only 
specular space harmonics. This is the familiar case of reflection and transmission 
from homogeneous layered media. 
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3 Numerical Solution for a Wedge- Corrugated Ab- 
sorber 


The above coupled wave method has been applied to the case of a wedge-corrugated 
absorber. The absorber is located in free space, and comprised of an iron-loaded 
epoxy material with approximate dielectric parameters e = 9 — j0.4 and /z = 1 — j 0.5 
at millimeter-wave frequencies [26]. The dimensions d and A are comparable to the 
free space electrical wavelength A = 2xfki. 


3.1 Discretization Criteria and Power Conservation 

The accuracy of the coupled wave solution depends on both the number of layers L 
and harmonics N. The effect of finite N is twofold: First, the dielectric interface is 
not modelled by a jump discontinuity in the material constants, but rather by a non- 
localized oscillatory transition characteristic of an incomplete Fourier reconstruction 
of a step functions. Second, sudden x— directed spatial variations in ~E and H are 
not described. To gain insight on the necessary value of N, we model the unknown 
field structure in a particular layer by a piecewise-continuous analytic profile. The 
profile (Fig. 5) is sinusoidal outside the absorber material with a period equal to the 
incident wavelength A. Inside the material the field decays exponentially with a decay 
constant determined by the material’s skin depth \/k%. 

A measure of the ability of a finite series of space harmonics to describe this 
profile can be defined by the ratio of energy in the lowest 2N + 1 Fourier coefficients 
to the total profile energy. For the model profile, the fractional power versus N is 
plotted in Fig. 6 for three layer positions: z x = 0,-d/2 and -d. When the composite 
field is primarily exponential in form (corresponding to z x near zero) a larger number 
of coefficients is required to achieve an acceptable fraction. For the above absorber 
material and with d = 7.5A and A = 3.75A (several electrical wavelengths), N > 
20 — 25 is required to obtain a fractional energy ratio of ~ 90 — 95%. This minimum 
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required value increases inversely with skin depth and wavelength. However, as N 
increases, the condition numbers of the propagator matrices increase, and numerical 
quantization noise begins to render the solution unstable. 

Electromagnetic power conservation requires that the incident z- directed power 
flux is equal to the sum of the z-directed power fluxes scattered into regions 1 and 3 
and the total power absorbed in region 2. The scattered fluxes are determined by 
integrating over a full period and representative width tv: 

fo Jo \ Jte{Ef(*. 2 = o) X 77r*(x, z = 0)} • Z dx dy 
ff / 0 A i Re(Ez (x, z = -d) x 7 7,-(x, z = -d)}- (-*) dx dy 
The power absorbed in the / th layer of region 2 is: 

r r / A Im. [ta{ £ ,(i, 2 ))ir a (i,2)l J + dli V iz 

i t|+1 Jo Jo 2 (50) 

Calculations based on Eqs. 48-49 show agreement with power conservation to within 
one part in ~ 10 6 . These calculations serve to detect the effects of numerical quanti- 
zation noise and illustrate the accuracy of the numerical implementation. The power 
conservation results are independent of the number of harmonics used in the analysis. 

Convergence of the total reflected power is monotonically increasing, and oc- 
curs near N fa 20 - 30. This is expected, since harmonic numbers at least this large 
are required for the dielectric interface model to appear discontinuous on the scale of 
field variation. For smaller N, the dielectric interface appears “fuzzy” on this scale, 
thus yielding a decrease in the reflected energy. The phenomenon is analogous to the 
decrease in reflectivity provided by a tapered matching section between transmission 

lines of different impedance. 

Due to the suitability of complex exponentials for describing both oscillatory 
and damped field variations in z, the required number of layers is not restrictive. 
Typically, discretization with Az w A/10 provides stable reflectivity values. 


(48) 

(49) 
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3.2 Graphical Visualization of Fields 

A graphical display of the fields around and within the wedge absorber has been 
developed to provide a qualitative verification of the accuracy of the coupled wave 
solution. Figs. 7a and b show electric field streamlines for TM incidence with #» = 30 , 
d = 7.5A and A = 3.75A. Discretization was performed using L = 19 and N = 25; no 
numerical quantization noise is apparent. The shadow region below and to the right 
of the wedge peaks can be clearly seen in Fig. 7a, consistent with geometric optics 
theory. Little specular reflection from the illuminated wedge face is apparent. This is 
expected since the angle of incidence relative to the face is close to Brewster s angle 
for TM waves 1 . 

A magnified view of the wedge tip region (Fig. 7b) reveals that the tangential- 
field boundary conditions are satisfied and that the field attenuation within the ab- 
sorber is consistent with the skin depth. Furthermore, there is evidence of surface 
polarization charge p, v near the wedge tip, where p tp = -t 0 n • [P \ - P%\- Here, 
P(x,z) = e 0 [e((x, z) - 1] £«(*,*) is the electric polarization vector for outside (1) 
and inside (2) the absorber surface and n is the outward-pointing surface normal. 
The p, p on the two opposing faces of a wedge forms a dipole surface charge wave that 
travels down the wedge. This dipole surface charge wave acts as a travelling- wave 
radiator, and is the principle cause of scattering from the wedge peaks. 


3.3 Comparisons with Integral Equation Formulation 

Results from the coupled wave formulation and an integral equation formulation using 

a periodic Green’s function (PGF) have been compared for consistency. The PGF 

formulation, detailed elsewhere [7], determines the fields in each region by numerically 

solving four integral equations, one each for the electric and magnetic vector potential 

functions on both sides of the wedge boundary. Because of the Floquet condition, 

l Since fji/rft is complex, there is strictly no angle of total transmission. However, the specular 
TM reflection coefficient exhibits a minimum at ~ 70° local incidence. 
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the integration is performed over one wedge period. The integrals involve unknown 
equivalent electric and magnetic currents defined along the wedge boundary, and 
known free-space Green’s functions for both sides of the wedge boundary. The method 
of moments is used to compute the unknown currents, and subsequently, the diffracted 
fields. 


The amplitudes of the zeroth and first-order reflected and transmitted mag- 
netic fields obtained using each of the two methods are shown for a non-absorbing 
periodic wedge surface (e 3 = 2.5, p 3 = 1.0) in Fig. 8. An incident TM plane wave 
with 6 X = 30° and A = 1.0A is used. In all cases, the illustrated amplitudes represent 
propagating waves. Both approaches agree well for depths up to 1.5A. However, the 
coupled wave technique has the advantage of being directly applicable to inhomoge- 
neous dielectric substrates. 


4 Laboratory Scattering Measurements 

Millimeter-wave bistatic scattering measurements of two wedge-corrugated absorber 
samples were performed to verify the coupled wave theory. Both absorbers consist 
of Emerson- Cumming CR-114 iron-impregnated epoxy cast to the dimensions in Ta- 
ble 6. The bistatic measurements were made using a focussed-lens scattering range. 
Deconvolution of the angular sensitivity pattern of the scattering range was performed 
using the CLEAN algorithm. 


4.1 Measurement System 

The scattering range consists of identical focussed-lens transmitting and receiving 
antennas interfaced to a Hewlett-Packard model 8510 network analyzer system 9. 
The HP 8510 system consists of a display /processor unit, two HP 8341 phase-locked 
microwave sources and an external S-parameter test set. The signal from the first 
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(transmitter) source is amplified and up-converted by a W-band mixer, then radiated. 
A directional coupler permits a portion of this signal to be downconverted and used 
as a reference for complex S 21 measurements. The received signal is similarly down- 
converted using the phase-locked signal from the second (receiver) source. The system 
provides a dynamic range of ~ 55 dB over the band 75-100 GHz. 

The two focussed-lens antennas consist of standard gain horns and rexolite 
lenses. Each are mounted on radial booms that rotate about a common axis lying 
in the plane of the sample’s surface. The configuration produces a focussed spot 
with a diameter of ~ 6 cm on the sample surface. The range permits both co- 
and cross-polarized bistatic scattering measurements for both TE- and TM-incident 
polarizations. 

Bistatic reflection data were measured for a range of 9{ from 20° to 70° in 10° in- 
crements. Due to a minumum 40° angular separation between each boom (for mechan- 
ical reasons), the corresponding range of 9 t was constrained to be 40° - 9 X < 9, < 70°. 
Although both magnitude and phase measurements were made, positioning errors 
rendered the phase information unusable. Bistatic data were also measured for a flat 
aluminum plate; this data was subsequently used for calibration and deconvolution 
purposes. 


4.2 CLEAN Algorithm 


The measured reflectivity amplitude S M is considered to be a convolution of the sam- 
ple’s angular reflectivity pattern S s and the angular response S p of the measurement 
system: 


SpJ6„4>,\6„4>.) = f’ jf* /! smff.iff, d#, 

° ° (51) 

where a,0 denote the transmitter and receiver polarizations (respectively), and are 
either TE or TM. Due to the relatively broad angular response of the measurement 
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system, the effects of S p must be considered. For identical transmitter and receiver 
polarizations, S p is obtained from the flat aluminum plate measurements. 2 Due to 
a lack of usable phase information, we have modelled the above observation process 
as a superposition of power contributions: 

dO'.d#, 

(52) 

In the above, we assume that contributions to the received power from all angles are 
incoherent. 




In order to compare the bistatic scattering pattern to the computed Floquet 
harmonic amplitudes, deconvolution of the probe function |S P | 2 was performed us- 
ing the CLEAN algorithm [21, 22). Developed for use in radioastronomical imaging, 
this iterative technique is well suited for reconstruction of point-like source functions 
from spatially low-pass filtered measurements. Since the angular reflectivity of the 
absorber samples consists of a discrete set impulses (each corresponding to a prop- 
agating Floquet harmonic), CLEAN is particularly useful. Further, CLEAN allows 
deconvolution of space- varying probe functions, as encountered here. 

The algorithm begins by initializing the cleaned array \S s \ 2 with zeroes, and 
the residual and probe arrays with |S M | 2 and |S P | 2 , respectively (Fig. 10). The peak 
in the current residual pattern is found, and a scaled probe pattern is subtracted from 
the residual at this location. The scaling factor is an empirical loop gain parameter 
selected to be in the range 0. 1-0.5. The remainder from this operation replaces the 
residual pattern, and a positive impulse of amplitude equal to the scaling factor 
is added to the CLEANed array at the location of the peak. The loop continues 
until all peaks in the residual pattern are smaller than the standard deviation of the 
measurement noise amplitude, which was 50 dB below the peak specular reflectivity 
from the aluminum plate. 

a Due to sensitivity limitations, our measurements were restricted to cases of like polarization. 
For differing transmitter and receiver polarizations, S p could be obtained using a quarter-wave plate 
in series with the flat reflector. 
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5 


Results and Applications 


5.1 Comparison of Measurements and Calculations 

The deconvolved reflectivities for the J and K samples were compared to coupled 
wave calculations of the power in the reflected Floquet harmonics. The comparisons 
were performed at a midband frequency of 89 GHz (A = 3.37 mm). The calculations 
used c = 9 - j0.4 and n = 1 - j0.5. The results of the comparisons are shown in 
Figs. 12 and 11 for TE and TM incidence, respectively. 

In most cases, the CLEANed data (|S 5 | 2 ) shows both specular (n = 0) and 
diffuse (n ^ 0) Floquet harmonics at approximately the expected angular locations 
(as determined by Eqs. 6 and 7). However, in some cases (due to measurement errors) 
CLEANing results in either additional (spurious) harmonics (Figs. 12b and 11a) or 
missing harmonics (Fig. lib). 

Referring to Eq. 45, the calculated reflectivity for the Floquet harmonic 
is: 

r„ = (53) 

with Tji H 0 = 1. As expected for surfaces with impedances tj 3 < ijj, the TE harmonic 
amplitudes are significantly larger than the corresponding TM amplitudes. Accord- 
ingly, TM is the desired polarization for best performance from wedge absorbers 
having e' > p'. 

The CLEANed reflectivities generally corroborate the coupled wave calcula- 
tions, although the measurements are usually higher by 3-10 dB. A large portion 
of this discrepancy can be ascribed to measurement accuracy. Although the short- 
term repeatability of the measurements was better than 0.5 dB, fluctuations of up to 
~ 3 dB were observed to be caused by mechanical misalignment. 
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5.2 Sensitivity to Assumed Dielectric Parameters 


We can determine the uncertainty in the calculated harmonic reflectivities using com- 
puted sensitivity functions: 


(*■-)• = + (&)'<■*? + (^) !<v)i + (w) 


wr 

(54) 


The 6c\ Sc", 6fi" are uncertainties the absorber material’s dielectric parameters, 
and are each estimated to be ~ 0.2. For our samples and incidence angles, the third 
term accounting for uncertainty in y! dominates in both TE and TM cases, 3 resulting 
in Sr nQ < 2 dB for any n. This uncertainty is hypothesized to account for some of 
the discrepancy between the computed and CLEANed reflectivity amplitudes. 


5.3 Wideband Absorber Design 


The design of wideband absorbing structures is based on minimizing the total reflec- 
tivity r 0 , defined as [27]: 


( 0 „ 4>i) = J- P* I * /2 53 sin 6 * de ' d< t>» ( 55 ) 

47T Jo Jo a 


where a,fi are either TE or TM and denote the observation direction. In the 
above, 4>i‘, 4 >>) is the surface bistatic scattering function for incident radiation 

of polarization a scattered into radiation of polarization fl. For periodic surfaces 7 ^ 
consists of a finite sum of delta functions, one for each propagating reflected harmonic. 


In terms of Eq. 4: 

cos 0 

7*»(*» *;*«*«) = 4jr L £(cos0.-cosM^.-WI^-«/H 2 -^^ ( 56 ) 

where the incident wave is of polarization a and amplitude E 0 - 1 , and e 0 is a unit 

vector in the direction of the electric field for /9-polarization. It follows from Eq. 55 

3 For the TM case, the coupling of the incident wave into the dielectric is more strongly influenced 
by /j' than (' since 77 is parallel to the edges. For the TE case, the opposite would be expected. 
However, since (' ss 9 and ft' « 1, the relative effect of changing c' by 0.2 is smaller than that of 
changing n' by 0.2. 
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that: 


(57) 


r„(»i,*) = E l*»l ! 


COS0, n 

COS 


For a perfectly reflecting surface, r Q = 1. 


Reflectivity calculations for a wedge absorber with normal TM incidence (0, = 
0° ), c = 9 — j'0.4 and fi = 1 — j0.5 and varying absorber dimensions have been 
computed using the coupled wave technique (Fig. 13). The number of harmonics and 
layers was 27 and 120, respectively. The results show generally decreasing reflectivity 
as the electrical depth of the structure increases, although the reflectivity tends to 
saturate for depths greater than ~ 1 — 2 A. This can be explained as follows: In- 
creasing depths increase the length of the inhomogeneous transition region between 
free space and the bulk absorber. This improves the match between these regions in 
the same manner provided by a tapered matching section between transmission lines 
of different impedance. The saturation is hypothesized to be the result of scattering 
by the tips of the absorber wedges. This mechanism appears to impose a lower limit 
on the reflectivity of a wedge periodic structure. 


In addition to the overall decreasing reflectivity, some periodic behavior can 
be seen, particularly for d/A from 0.25 to 1.0. This behavior can be explained by 
noting that sharp increases of reflectivity tend to occur at integral values of A/A, 
corresponding precisely to the onset of additional propagating Floquet harmonics. 
Such harmonics carry additional power away from the interface. In the limit of a flat 
surface, the reflectivity converges to the specular value (-6.2 dB), as required. 


The calculations suggest that wedge-type absorbers with total TM reflectivities 
less than -30 dB (0.001) and decade (or greater) bandwidths can be obtained by 
choosing d/A > 2. Although these calculations assumed specific dielectric values, the 
reflectivities for absorbers of slightly different dielectric constitution can be estimated 
using the sensitivity parameters in Table 6, where: 


r„(e + 6£,p + «f0»r„(e,|i) + + §£& " + V + (dB) (58) 


dfi' 
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The listed sensitivitiy parameters were obtained numerically by computing divided 
differences. In the above, d' and fi" are positive for lossy media. 


6 Summary 


Linear N-port analysis along using coupled wave theory provides an efficient numerical 
method for computing the fields both within and scattered from a lossy periodic 
structure. Reflectivities computed using this method compare favorably to measured 
laboratory data after deconvolution using the CLEAN algorithm. Discrepancies in 
the comparison are ascribed primarily to measurement uncertainty. The coupled 
wave reflectivities compare very well with calculations based on an integral equation 
formulation using a periodic Green’s function. A distinct advantage of the coupled 
wave method over Green’s function methods is its ready applicability to arbitrary 
periodic grating structures and otherwise inhomogeneous periodic dielectric materials. 
This includes structures consisting of several corrugated layers of different material 
and periodic dielectric structures containing voids. 

In principle, the coupled wave method is applicable for arbitrary dielectric 
material parameters. However, for materials with skin depths that are much smaller 
than the wavelength (e.g., good conductors), a large number of harmonics are required 
to describe the field variations. In this case, the the coupled wave method requires 
the inversion of large ill-conditioned matrices, resulting in numerical noise in the 
solution. This limitation can be surmounted by higher arithmetic precision, but at 
increased computational expense. Alternate noise mitigation methods include least- 
squares boundary matching techniques (e.g., [28]) and ad-hoc elimination of decayed 
harmonics from layers in which they are known a-priori to be negligible [20]. 

Using the coupled wave method the reflectivity of representative wedge-corrugated 
absorbers has been computed. By proper selection of the pitch and depth the total 
TM reflectivity at normal incidence can be reduced to lower than -30 dB over a decade 
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in bandwidth. It is the total reflectivity, rather than the retroreflectivity, that is im- 
portant in the design of microwave calibration targets. Additional requirements are 
that the temperature distribution within the target be precisely known; the problem 
of the temperature distribution within a microwave calibration target and its effect 
on the thermal emission are being addressed separately. 

Some straightforward extensions of the coupled wave theory follow from the 
above formulation. These includes the obvious cases for which the grating vector is 
not in the plane of incidence (<f> ^ 0) and for which the surface is periodic in two 
dimensions. Other innovative applications include scattering from dielectric extru- 
sions using a coupled cylindrical wave formulation and scattering from an arbitrary 
dielectric object using a coupled spherical wave formulation [29]). 
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Table 1: Dimensions of wedge-corrugated absorber samples used in the laboratory 
bistatic measurements. 


Sample 

A 

d 

J 

3.2 

6.5 

K 

12.5 

25.4 


Table 2: Sensitivity of total TM reflectivity with respect to dielectric parameters 
(dB/unit) for normal incidence. The mean dielectric parameters are c = 9 - j0.4 and 
H = 1 - j0.5. 


Quantity 

A = 
d/A = 1 

0.5A 
d/A = 4 



II 

< 

2A 

d/A = 4 

ttm (dB) 

-35.1 

-37.7 

-27.3 

-40.9 

-28.6 

-33.5 

drjM 

dt f 

0.90 

0.25 

0.68 

-0.38 

0.45 

0.45 

drjM 

de" 

-4.01 

0.14 

0.09 

-0.30 

-0.72 

0.29 

drjM 

dfi' 

-68.6 

-3.78 

-6.09 

-10.57 

-11.4 

-2.70 

drjM 

dn " 

-11.7 

4.66 

-6.88 

2.77 

-4.61 

2.69 
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Figure 1: Periodic surface scattering geometry for an obliquely incident, arbitrarily 
polarized plane wave. 

Figure 2: Schematic diagram illustrating the implications of the Floquet condition 
and phase matching on propagating and evan ascent plane waves in the upper and 
lower half spaces. 

Figure 3: Discretization of the periodic region (2) into N layers for the coupled-wave 
analysis. Each layer is modelled as a linear multi-port network with an associated 
transmission matrix. The scattering problem is solved by cascading the individual 
layers. 

Figure 4: Schematic diagram of a single grating layer, illustrating incident and scat- 
tered waves. The scattering behavior of each grating layer is described using a trans- 
mission matrix by considering the layer to be a linear multiport network. 

Figure 5: Piecewise-continuous analytic profile used to model the field structure in a 
period layer. 

Figure 6: Fractional power in the incomplete Fourier series representation of the 
piecewise-continuous analytic profile. 

Figure 7: Electric field streamlines for a TM plane wave obliquely incident on a 
wedge absorber: (a) full view, showing geometrical shadow (below dashed lines) and 
(b) magnified view of a wedge tip, showing surface polarization charge density. The 
parameters are 0, — 30°, d = 7.5A, A = 3.75A, c = 9 — jO.4 and fi = 1 — jO.5. 

Figure 8: Comparison of transmitted and reflected harmonic amplitudes (arbitrary 
units) for each of two methods: coupled wave and periodic Green’s function. The 
magnetic field magnitudes for the indicated propagating harmonics are plotted. 

Figure 9: Schematic diagram of the focussed-lens millimeter-wave scattering range. 

Figure 10: Functional block diagram of the CLEAN algorithm. 

Figure 11: Comparisons of measured data and coupled- wave calculations for wedge 
absorbers with TM incidence and 0, = 60° : (a) J sample, and (b) K sample. 

Figure 12: Comparisons of measured data and coupled- wave calculations for wedge 
absorbers with TE incidence and 0, = 60° : (a) J sample, and (b) K sample. 
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Figure 13: Wideband TM reflectivity vs. period A for wedge absorbers of varying 
depth d. The calculations used the coupled wave method. 
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